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It is shown that the universally unstable kinetic drift wave in an electron-ion plasma can very 
effectively be suppressed by adding an extra flowing ion (or plasma) population. The effect of the 
flow of the added ions is essential, their response is of the type ( v p h — f/o) exp[— (v p h — vgo) 2 ], where 
Vfo is the flow speed and v p h phase speed parallel to the magnetic field vector. The damping is 
strong and it is mainly due to this ion exponential term, and this remains so for vgo < v p h- 

PACS numbers: 52.35.Kt, 52.30.Gz, 52.35.Qz 


I. INTRODUCTION 

Drift wave is called universally growing mode due to 
the fact that it is unstable in both fluid and kinetic de¬ 
scriptions, and in collisional and collisionless plasmas. 
The wave is self-excited and it grows due to the free en¬ 
ergy in plasma inhomogeneity and this remains so even 
in plasmas with hot ions. This fact was used in our re¬ 
cent papers where a new paradigm was put forward for 
the heating of the solar corona [2 |2j and solar wind [3j 
by drift waves, based on stochastic heating mechanism 
from Refs. HIS]. In practical situations in lab plasmas, 
in order to study the mode in a controlled situation, the 
wave is on purpose excited and driven by an electron cur¬ 
rent [B] or by a shear flow [Tj. More recent experimental 
studies of drift instabilities are available in Refs. 00. 

This is a dangerous mode in any plasma environment 
pTOj-pS], and various effects have been studied in the past 
in order to stabilize it. One of them is the magnetic 
shear which in simple slab geometry introduces a layer, 
in the direction of the shear gradient, at which the mode 
is stabilized by resonant ions HU HU. So although the 
stabilization is kinetic by nature it is routinely described 
as an effect of the plasma geometry. 

But in more realistic laboratory situations, the same 
geometry which implies the stabilization by the mag¬ 
netic shear in fact includes some additional features, like 
toroidal mode coupling, which may completely cancel the 
magnetic shear stabilization [18]. Much more on these 
phenomena may be found in our earlier works (19H21) . 

Yet another way of the drift wave stabilization is by 
cold electrons added to the plasma. In the present work 
we show that this can also be done by adding flowing 
ions or plasma , which need not be cold at all. 


II. THE MODEL AND DERIVATIONS 


dex / further in the text). We allow for the presence of 
electrons as well in the /-species in order to avoid the is¬ 
sue of excess charge in case that ions alone are added, i.e., 
n e /o = n.ifo = const., although stabilization is mainly by 
the /-ions. The equilibrium density gradient is in the 
x-direction and the wave vector k lies in the y , z-plane. 

Electrostatic perturbations are assumed propagating 
nearly perpendicular to the magnetic field ~ exp (—iuit + 
*kj_r + ik z z). The perturbed densities are: 
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The geometry assumed in the derivation is such that 
the background magnetic field (B 0 ) is oriented along the 
z-axis. We assume a static (denoted by index s) inhomo¬ 
geneous and quasineutral electron-ion plasma n es o(x) = 
niso(x), penetrated by a homogeneous plasma stream (in- 
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a denotes the species, q a is their charge, n a o is the equi¬ 
librium density, L nsa is the inhomogeneity scale length 
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of static component and u; = ui — k z v a o is the Doppler 
shifted frequency due to the streaming velocity v a o- 

We are considering the case of low plasma beta j3 a = 
2fj 1 Qn a oT a /B q « 1 due to which the magnetic field gra¬ 
dient is ignored following the relation L nsa /Ls S a ~ (3a 
j3], where Lssa is the scale length of magnetic field in¬ 
homogeneity. The parallel integration gives rise to the 
plasma dispersion function W^(£„( Sl /) a ) with the argu¬ 
ment £, n (sj) a , where the perpendicular integration yields 
the modified Bessel function in the term An( 6 ( S) /) a ) = 
e -b (s./)<» I n (br at f) a ) with the argument br s ,f)a‘ For the 
static component r) SQ o = 0, u = ui, and w* e , u>* si ^ 0 , 
while for streaming particles Vf a o ^ 0 and u; = w — 
k z Vf a o, oj*f a = 0. The drift frequencies for electrons and 
ions are related as u>* si = — (T si /T se ) ui* e where w* e > 0. 

The dispersion relation for the electrostatic drift waves 
is obtained from linearized Poisson’s equation 

£ok 4*1 — {'('H'esl T Tlefl Hisl l)- 


The Larmor radii of electrons in both plasmas are very 
small as compared to the ions, which allows for the ex¬ 
pansion of the modified Bessel function for small argu¬ 
ment as A n ( 6 ( s j) e ) = [6( s j) e /2] n /n\. It is easy to see that 
only n — 0 terms survive in the limit of a negligible value 
of the argument, i.e., for 6( s j) e — ► 0, Ao(&( s ,/)e) = 1- We 
shall also separate the n = 0 terms in the ions contribu¬ 
tion. 

Using the identity A n (x) = A_ n (x) and the expansion 
of the plasma dispersion function for n / 0 terms in limit 
of the large argument, and assuming the realistic low 
frequency case for both the components, i.e., £, n fi, £,nsi ^ 
1 and ui < f 2 ( Sl /)j, respectively, one can easily prove that 
the n/ 0 terms vanish from the last terms of Eqs. 00 
and we get the dispersion relation 
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In order to calculate the growth rate of the drift wave, 
we separate the real and the imaginary parts of Eq. ([5]) 
and the growth rate becomes 
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The real dispersion relation may be obtained by taking 
e r = 0, i.e., 

3 + Re [A s + Af] / (k 2 \ 2 Dse ) = 0. 

The wave behavior will be discussed in two different fre¬ 
quency limits. 

A. Specific frequency limits 

a. In what follows the electrons and ions in the 
static component satisfy the following frequency limits 
k z Vtsi € w € k z Vtse, while for the streaming species we 
have 


k z v t fi < |w - k z v f0 \ < k z v tfe . (9) 


The growth rate for the drift wave becomes 


7i = —Ci (51 + g 2 +£ 3 ) = 


/tT ui 2 n es o/n is0 

V 2 w* e A 0 (b si ) 


X 




+Aq (b S i) 


'^'iso'^’se M r M 
^eso'-^’si kz^tsi 


exp 


^ r V 
2 k 2 vli)_ 


HifO^se ^r k z VfQ 

TlesO Tfe k z Vtfi 


A o(bfi) 


Tfe 

Tfi 


exp 


- k z Vf 0 ) 2 \ 

2k z V tfi J 



exp 


(U> r - k z Vf 0 ) 

2k 2 z v 2 fe 



( 10 ) 


Here, the meaning of the terms < 71 , < 72 , <73 is obvious, and 
they describe contribution of static electrons and ions, 
and flowing electrons and ions, respectively. The real 
part of the frequency may be written as 
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Here, Ps(s,f) /Hq c^fs,/) The 

static ions term g 2 in Eq. ( flO] ) cause damping regard¬ 
less of parameters, while the electrons term yields the 
usual kinetic instability provided that necessary condi¬ 
tion uv < oj* e is satisfied. 

As for the contribution of the flowing plasma, the ( 73 - 
term , it turns out that the universally growing mode can 
completely be stabilized and this will be demonstrated 
below using some parameters that may be applicable to 
the laboratory plasma conditions. 

We choose parameters which will show that damping 
by /-plasma is essentially due to their flow. We take 
So = 2 T, L nsi = L nse = L n = 0.1 m, T si = T se = T s = 
10 5 K, riiso = n es0 = n s0 = 10 19 m^ 3 , A_l = X y = 3 mm, 
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FIG. 1: The imaginary part of the drift wave frequency (101 
normalized to u} r in terms of the flowing plasma density, for 
the same temperature of both plasmas. 



FIG. 3: The flowing plasma term g$ from Eq. (101 used in 


Figs. [T] [2] in terms of the /-plasma speed Vfo . The peaks are 
the /-ion exponential term contribution. 



FIG. 2: The imaginary part of the drift wave frequency (101 
normalized to oj r in terms of the /-plasma speed, for the same 
temperature of both plasmas. 



FIG. 4: The growth rate/damping (101 for a cooler /-plasma 
Tf/T a = 1/5 in terms of the speed Vf 0 . 


and take k z /k y = 0.0002. For such parameters the drift 
wave is unstable, ui r = 74659 Hz, ji/oj r = 0.034 in spite 
of so hot s-ions. When /-plasma particles are added, and 
with the same temperature Tf z = Tf e = T s , there is very 
little change in 7 1 and the wave remains growing even if 
n/o is strongly increased, and this remains so as long as 
/-particles do not flow (see full line in Fig. [l]). 

But if the /-plasma is flowing, the growth/damping 
is changed. In this case the initial instability caused by 
the electron term g\ is first increased for small nfo, but 
for larger /-plasma density the mode is heavily damped. 
The effect of the flow of the added plasma is thus es¬ 
sential, it has a profound effect on the drift wave. See 
more details in Fig. [2] which shows that it is possible to 
find particular speed values for which the mode is most 
effectively damped. Note that here v/o < Vph = to r /k z = 
148530,137104,127311 km/s for n f o/n s0 = 0.2, 0.3, 0.4. 

The effect of the flow may be understood from Fig. [3j 
where the /-plasma term 53 is presented in terms of Vfo 


for the three values rifo ■ The essential part is the ion 
term, which is of the shape (v p h — Vfo) exp[— [v p h — /o) 2 ] 7 
so that normally destabilizing first part v p h — Vfo is coun¬ 
teracted by the ion exponential part, and the /-ion part 
in g 3 goes to zero for large v instead of linearly increas¬ 
ing the growth rate indefinitely due to ( v p h — u/o) term 
alone (for Vfo > v p h)- The electrons have a minor role 
and contribute only in the range Vfo > v p h when the flow 
destabilizes the mode (not presented here). The lines are 
made broken to mark regions where k z Vtfi <C |w — k z Vfo\ 
is violated and this analytical model should not be used, 
hence the given speed limit here and in Figs. ii 

To check the effects of the /-plasma temperature, in 
Fig .[4] we give the imaginary part of frequency in terms of 
Vfo for the two densities n/o/^so and for Tj /T s = 0.2. It 
is seen that for n/o/^so = 0.1 the mode is immediately 
damped even for Vfo = 0 as soon as the /-plasma is 
added, and there is a strong damping for larger u/o- The 
dependence of 71 on Tj is complicated, but Fig. [4] may 
partly be understood from Fig. [ 5 ] where we set n/o/n s o = 
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FIG. 5: Flowing plasma term g 3 in the imaginary part of FIG. 6: Drift wave stabilized by flowing ions, 

frequency (| 10[) in terms of /-plasma temperature. 


0 .1; it is seen that in the given temperature range, around 
Tf/T s = 0.2, the g 3 term in (10) is positive and it causes 
strong damping, but this is not always so. 
b. In the frequency range 


\u - k z v f0 \ « k z v t fi, 


( 12 ) 


the /-ion response is nearly Boltzmannian while electron 
contribution is negligible (see further in the text), and 
with similar approximations as above we have: 


our Vtsei here the right-hand side can clearly be much be¬ 
low unity, so the ion flow in this regime can be far more 
efficient in exciting the drift mode. 

The omitted electron terms make only minor changes 
in Eqs. (13 14): a 3 term is multiplied by a small factor 
l + (m e /mi)^ r? (Ty i /T/ e ) 3 / 2 , and the last term in denom¬ 
inator of Eq. (14) is multiplied by a term 1 + Tfi/Tf e . 


B. Flowing ions case 
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Here, n es3 = n es0 = n s0 , n ef0 = n if0 = n f0 ■ The /- 
ion term a 3 causes a strong damping when vjq is small 
and this can be checked for the same parameters as be¬ 
fore. However, for Vfo > v p h the wave is destabilized 
and this can easily be more efficient than in the case 
of electron-current driven mode bei. Indeed, in the 
usual electron-ion plasmas, the latter implies an addi¬ 
tional electron current term oq = uo/v Te in the growth 
rate but this can easily be smaller than the exist¬ 
ing a 3 term. For Vfo > v p h we have that a 3 > u,o/v Te 
if Vf 0 /u 0 > ( Vtfi/v Te )(Tfi/T se )(n s o/nfo ). Taking v Te as 


We checked the case of adding flowing ions only, in 
the range k z v t fi \lu — k z Vfo\, assuming that plasma 
adjusts in such a way that global quasineutrality is pre¬ 
served n es 0 = ni S 0 + n-ifo- This is completely equivalent 
to Ref. [22. where the stabilization is discussed by an ad¬ 
ditional cold electron population. In Eqs. (10 11) vanish 
the /-electron term, and the factor 1, respectively. The 
result is presented in Fig. [6] for several densities of the 
flowing ions and the result is similar to Fig. [2] The fre¬ 
quency is oj r = 70926,67193,59723 Hz for nf i0 /n se0 = 
0.05,0.1,0.2. Here we keep Tfi = T si = T se = 10 5 K, 
and other parameters are the same as before. The wave 
behavior is very similar to the previous plasma flow case. 

Here, the /-ion flow in principle implies a current that 
might cause a sheared magnetic field component B s = 
g,oenfoVfoL s , where L s is the characteristic shear length, 
which is known to stabilize the drift wave itself da Eg. 
However, for parameters used in the text the sheared 
component is negligible; at the perpendicular distance 
L s = L n it remains below O.OOlHo- At shorter distances 
it is even smaller and can be neglected. 


III. SUMMARY 

In conclusion, this work provides some clear recipes 
for damping of the drift wave which is usually believed 
to be universally unstable. The stabilization is expected 
to work for various modes from the drift wave spectrum 
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and it can be used as an alternative for some other mech¬ 
anisms proposed in the past [24H26j . 

The model presented here also has an obvious advan¬ 
tage with respect to so called stabilization by cold elec¬ 
trons (having some temperature T c ) [22] because the lat¬ 
ter disregards simultaneous cold electron collisions with 
other species (which is proportional to l/T^ 2 , so the 
cooler electrons the more collisions). Hence, these colli¬ 
sions can be frequent even if the plasma is fairly collision¬ 
less regarding its usual components (ions and hot elec¬ 
trons). On the other hand, thermalization of such cold 


electrons is instant and its characteristic time is the same 
as their velocity relaxation, and they are thus totally inef¬ 
ficient in stabilizing the drift mode. So the ion (plasma) 
flow stabilization presented here is clearly a far better 
alternative. 
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